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^ . Abstract 

OO 

ff^ ' Perturbation theory for non-Abelian gauge theories at finite temperature is plagued 

y^ . by infrared divergences caused by magnetic soft modes ~ g^T, which correspond to the 

fields of a 3d Yang-Mills theory. We revisit a gauge invariant resummation scheme to solve 

this problem by self-consistent mass generation using an auxiliary scalar field, improving 



cn 

fvj . over previous attempts in two respects. First, we generalise earlier SU(2) treatments 

to SU(Y). Second, wc obtain a gauge independent two- loop gap equation, correcting 
an error in the literature. The resulting two-loop approximation to the magnetic mass 

ij , represents a ^ 15% correction to the leading one-loop value, indicating a reasonable 

%^ ' convergence of the resummation. 



1 Introduction 

From heavy ion collision experiments to astroparticle physics and cosmology, there is a need 
for theoretical predictions from finite temperature gauge theories within and beyond the stan- 
dard model. Since Monte Carlo simulations of lattice gauge theories do not work for finite 
baryon densities or dynamical problems involving real time, analytical approaches are war- 
ranted as well. Perturbation theory for equilibrium quantities of non-Abelian gauge theories 
at finite temperature T features three scales: a hard scale ttT of the non-zero Matsubara 
modes related to the compactified Euclidean time direction, the soft and ultra-soft scales 
gT, g^T, with gauge coupling g, which are associated with the screening of colour-electric and 
colour-magnetic gauge fields, respectively. Since the dimensionless expansion parameter for 
the latter features a mass (or momentum) scale in the denominator, ~ g^T/E{p), prohibitive 
infrared divergences occur in the perturbative series when bare, massless propagators are 
used. This is known as the Linde problem [HIS]- Since the dynamically generated magnetic 
screening mass is itself m ~ g^T, the problem cannot be cured to any finite order in ordinary 
perturbation theory. On the other hand, evidence from gauge fixed lattice simulations (see [3] 
and references therein) as well as gauge invariant simulations of field strength correlators [1] is 
consistent with an effectively massive gluon propagator. Moreover, a dynamically generated 
gauge boson mass plays a prominent role in a Hamiltonian description of three-dimensional 
Yang-Mills theory [5l |6] and is of renewed interest in the context of a Higgs-less electroweak 
gauge sector, see e.g. [7]. 

In this work we systematise a resummation method for 3d Yang-Mills theory, viz. the 
magnetic sector of finite T gauge theories, which has been proposed some time ago [8]. 
The general idea is to screen the infrared divergences by adding a gauge invariant mass 
term, which gets subtracted again at higher order, resulting in a resummation of the loop 
expansion. The mass term is not unique and several possibilities have been tried at one-loop 
[El [T0| [TT] . Under the name of screened perturbation theory similar techniques have been 
applied to scalar theories |121 [H] and to the colour-electric sector ~ gT of gauge theories in 
a dimensionally reduced setting |14| as well as in four dimensions |15) . There, a screening 
mass is generated in ordinary perturbation theory and the resummation of its corrections 
is merely used as a means to improve the convergence. By contrast, the magnetic mass 
~ g^T is entirely non-perturbative. It has to be generated by an infinite resummation and 
evaluated self-consistently by a gap equation. Since the gauge coupling drops out of the 
effective expansion parameter, there is no parameter to tune and hence no regime where 
the resummation is guaranteed to work. Its convergence properties can be judged only after 
explicit calculations. 

In this paper we present a systematic derivation of a gauge invariant resummation scheme 
using auxiliary fields based on the non-linear sigma model. In particular, we discuss the BRS 
invariance of the resummation scheme which ensures that the symmetries of the theory are 
maintained throughout a modified perturbative treatment. We compute the gap equation for 



the magnetic screening mass through two loops in general R^ gauges and duly find it to be 
gauge parameter independent, thus correcting an error in an earlier two- loop investigation 
|16j . The two-loop result amounts to a ~ 15% correction to the leading one-loop result, thus 
pointing at a reasonable convergence of the resummation scheme. 

2 Resummation based on a non-linear sigma model 



Let us now focus on Euclidean Yang-Mills theory. 

1 

where we use a matrix notation. 



^(^) = -TT72Tri^'-, (1) 



F^, = [Df,, D,] , D^ = Id,, - A^{x) , A^{x) = igT'^Alix) , (2) 

with hermitean SU(A^) generators T", a = 1 . . . A^^— 1. We are interested in three dimensions, 
fj.,1' = 1 ... 3, where the coupling constant g^ carries dimension of mass and is to be identified 
with g^ = g\^-, if the action is viewed as the magnetic part of hot Yang-Mills theory. 

The general idea of a resummation is to sum up higher order contributions (infinitely many 
in our case) into a given order of a perturbative expansion. In order to avoid double counting, 
these contributions must then be left out at the higher order where they naturally occur, such 
that the perturbative scheme gets reorganised in some systematic way. (For a discussion of 
various schemes used in the context of thermal field theory, see |17j). This can be formalised 
by rewriting the Lagrangian serving for the perturbative expansion as |10] 

1 r 



/:eff ^ 



C{^iX) + A£(V^X) - ^A£(\/ZX) , (3) 



where X generically stands for the fields. The modification A£ contains fields of the original 
Lagrangian and possibly auxiliary fields. In particular, if A£ is chosen to represent a mass 
term for the gluon, this will regulate the infrared divergences. The counting parameter £ in 
which one expands is to be set to £ = 1 at the end of a calculation, for which the Lagrangian is 
identical to the original one. However, in a perturbative evaluation to finite order the results 
will differ from the unresummed ones, the original theory being recovered exactly only at 
asymptotically high orders. Whether low order calculations are a good approximation to the 
full answer has to be judged from the apparent convergence of the series and may depend on 
the observable, just as in ordinary perturbation theory. 

A valid resummation scheme has to maintain the symmetries of our gauge theory at every 
order, for general L Clearly, this leaves many conceivable choices for A£, several of which 
have been tried at one-loop level in the literature [U El [ini HI]. An optimal choice would be 
based on convergence properties in higher orders. Here we work with a gauged non-linear 
sigma model, coupling a field <I> G U(A^) as 

2 

A£(y4, vr) = ^Tr [(Z)„$)ti?^$] , ^^^ = 1 , $(x) = e'^(^) , nix) = i^T^TT^ix) , (4) 
g m 



where the 7r"(x) G IR are auxihary would-be Goldstone boson fieldqll with the same mass 
dimension as the gauge fields A"^(x). Under gauge transformations, $' = U^ and (D^^)' = 
U{D^^), with the unitary matrix U = e^^' and corresponding real coefficients A"(x), with 
A(x) = iT°'A°'{x). Thus, AC provides a mass term for the gauge fields at tree- level while 
maintaining gauge invariance. 

3 Gauge fixing and BRS- invariance 

In order to do perturbative calculations, a gauge needs to be fixed. It is well known that in 
Higgs and sigma models the standard covariant gauges lead to non-diagonal or mixing terms 
in scalar and gauge fields, ~ (S^vr)^^. This can be avoided by choosing i?g-gauges. In the 
case of a resummed calculation, however, additional choices have to be made. 

We can either take the point of view that our starting point is the resummed theory before 
gauge fixing as in Eq. ^, and then add gauge fixing and ghost terms to that expression, 

J^eSA = j{c{V£A) + {l-l)[A/:iVeA,ViTr) 

+£gi^A{V£A,ViTT) + CFPAiV£A,ViTT,Vec)j^ , (5) 

C^fMATT) = -^Tr!^[id,A,)-Cm'Tr[i^-^^)T^)T^yy (6) 

£fp,a(^, 7t, c) = \tv [2{d^c) {{d^c) - [A^, c]) + ^m^c (<^^c + c$) } , (7) 

with ghost fields c = igT"'c"' and anti-ghost fields c = —igT"'(?'. We refer to this gauge fixing 
procedure as A. For £ = 1 we obtain Yang-Mills theory without gauge fixing, as in Eq. ([3]). 
This invokes the following Feynman rules for the counter term two-point vertices, 

-♦- : ^tuA^') = (^ + mX^) M'^" , (8) 

_»_ : Tf{^^) = {p' + im^y5''\ (9) 

■■►■•■►■■ : Tf{c^) = {p^ + im^yS'^K (10) 

Alternatively, we may consider Yang- Mills theory in a covariant gauge and then resum the 
gauge fixed theory, a strategy which we refer to as procedure B, 

Aff,B = j[c{^rtA) + {l-^)AC{ViA,Vi^) 

£gf,B(Avr) = -^T,\^({d,A,)-{l-l)im'TT({^-^^)T-)T-)'Y (1^) 
£FP,B(^,vr,c) = ^Tr {2(a^c)((V)-[^^,c]) + (l-^)Cm2c($tc + c$)} . (13) 



^Note that it is the vr and not the <I> who get rescaled as tt — ?> yl-K for the purpose of resummation, Eq. (|3|). 



The corresponding counter term two-point vertices read 

-w#-- : r»t,B(^') = [mH^^W , (14) 

r|''(7r2) = (/ + ^(2 - i)my6'''' , (15) 

Tf{c^) = (em^)^"^ . (16) 



The two formulations feature non-trivial differences. Note that Eq. (J15p contributes to both, 
order ~ i^ and ~ i'^. Let us remark here that the gauge fixing and corresponding counter 
terms used in [8l [16] work only for the gluon pole mass to leading and next-to-leading order, 
but require generalisation (as above) for higher orders and other observables. 

In order to provide gauge invariant results for physical observables order by order in pertur- 
bation theory, it is necessary and sufficient that the gauge fixed Lagrangian Cgs is invariant 
under BRS-transformations |18[ [TU] . Note that in the resummed theory with general £ all 
fields get rescaled by yZ, and so do the gauge transformations. A BRS-transformation now 
corresponds to the special choice A = coylc. The variations of fields and Faddeev-Popov 
ghosts under infinitesimal BRS-transformations are 



SbA^ = 


= uj{d^c) + ujVi{cA^-A^c) , 


Sb^ = 


= LoVic^ 


^ 




n=0 "• j=0 ^-^^ 




^ Loc (vrc — cvrj H (vrvrc 



- 2TTCTT + cvrvr) + ©(Tr^,^^) , 

5bc = cjJvIcc , 

6bc = -j(^{d^A^)-{l-i)Cm''Tr(^{e^'^-e~^'^)T-)TyVe), (17) 

where Bn are the Bernoulli numbers. The above transformation refers to CeS,B] for formula- 
tion A, the factor (1 — i) in 6bc is absent. 

As a non-trivial check we have performed our calculations in both ways, obtaining identical 
gauge invariant results for both setups. We present our results according to setup B, as it is 
closer to the standard perturbative treatment. 

4 Relation to SU(2) calculations 

Here, we connect our general SU(A'^) parametrisation of the scalar field, Eq. (j4j), to the special 
case of SU(2) treated in P[l6]. Using Tr (T'^) = 0, Tr (T^'T^) = ij'^* and Tr 1 = iV, one gets 
for the product of two generators the standard expression 

r"r* = — i + Tr({T",T''}r")r^ + Tr([T'',r'']r^)T'= 

= ^1 + 2^'*^'^^'' +-i^^'=T^ (18) 



Note that for the special case of SU(2), where T" ~ a'^ and the Pauh matrices anticommute 
as {(t",(T } = 2(5 1, the totally symmetric structure constants vanish, d = 0. Hence, in 
SU(2) the product of two of our scalar fields is diagonal. 



\2„a ft 1 /^" T, I mhcrrc\ SU(2) /„• /r,_N2„a„a 



7r7r = (i5/m)'7rV-( -— l+(i'"'^rM ^' (ig/2m)'^''^n , (19) 

and therefore the field $ can be expressed as 

$ = e^ = cos- + isin- ^^^^ al+ivf'^r'^ (20) 

i i 

with a = cos iX^^ ~ 1 + o(^2) ^ (21) 

2m 



and ^° = , - s^in ^ ~ JL vr" + O TT^) . 

^/tt'^tt'^ 2m m 



(22) 



Hence, for N = 2 our model can be recast into the form of the model considered in [8]. For 
general N, however, the coefficients on the right-hand side of $ = -j^ 1 + 2Tr {^T"") T"" can 
not be expressed in a closed form in terms of the real fields vr". The traceless part of $ is 
part of our construction for the gauge fixing term, see Eqs. ([6]) and (I12p . which represents a 
non-trivial generalisation of R^ gauges to SU(A^). 

5 Pole mass from a gauge invariant gap equation 

Having designed a general gauge invariant resummation scheme for 3d Yang-Mills theory, we 
now apply it to a calculation of the gluon self-energy. Its transverse and longitudinal parts 
Hrp/L are defined as 

Kip) = '^"''{(v-^)nT(/) + ^H^(p^)} . (23) 

The self-energy itself generally is a gauge dependent quantity. However, the pole of the 
transverse part of a gauge boson propagator, Dti is known to be gauge invariant order 
by order in perturbation theory [201 I21j . and we may employ our resummation scheme to 
evaluate it. (The longitudinal degrees of freedom with a gauge dependent pole correspond 
to unphysical would-be Goldstone bosons as usual in Higgs-like theories, and can be gauged 
away in unitary gauge.) Without resummation the pole of the bare transverse propagator is 
at p^ = 0, whereas in the resummed theory it gets shifted to p^ = —rr?. Identifying m with 
the pole of the full propagator, we require that the pole stays at i? = —m? to any loop order. 
Taylor expanding the self-energy about the pole, the transverse propagator reads 

D^ = I P^=-™^+V ^-"t(--) .24) 

p2 + ^2_nT(p2) - i!^;r(p^l) +v + 0((V)2)' 



where n;^(— tti^) = dp2llT^{p'^)\p2^_^2 . Near the pole it then corresponds to a massive prop- 
agator with residue Z{m?), 

^T«4x^' ^("^') = r-fr7 — 2T' (25) 

pz _|_ j^z 1 — ilrp(— rri^j 

provided the first term in the denominator of Eq. (|24|) vanishes. This leads to the gap equation 

! nT(-m^) 
°= i-n-^(-m^) - ('^) 

Introducing the ^-expansion of the self-energy, IIt = Z^n>i ^"^rp" (p^), we now expand the 
gap equation to the desired order in £ and evaluate it after setting i = 1. Note that, to every 
order in £, the gap equation receives different kinds of self-energy contributions contributing 
to order £". In order to list these separately, we introduce nJp" to denote the sum of 

diagrams with (n — k) loops and k counter term insertions, 

n 

nf)=J^n5?"'^)'\ (27) 



The corresponding gap equation combines gauge dependent contributions of different self- 
energies nJp' into a gauge invariant quantity. 

We perform our calculations using dimensional regularisation, working in d dimensions and 
with generic gauge fixing parameter ^. Details and intermediate results are relegated to the 
appendix, from which we collect the results in the following sections. 

5.1 One-loop gap equation 

To leading order i^, the gap equation Eq. (f26]l is simply 

= <)(-m2) = n^^)'°(-m2) + uP'\-m') . (28) 

The five diagrams shown in Fig. [1] have been computed in a general i?^-gauge in [8] and lead 
to a .^-independent gap equation when evaluated on the pole. From appendix[Al we reproduce 
these results for SU(A^) as 

^niO/ 9n d=3-2e q^Nm /3 63 \ ,^, , , ^ 

nW'°(-m2) « ^__(^___in3J+0(6), (29) 

Iif^\p') = mK (30) 

Solving the quadratic one-loop gap equation Eq. ()28p then yields the well-known solutions 

mi_ioop = or mi_ioop= (^^ln3-^j ^ = 0.142276 5^ TV. (31) 



Figure 1: The self-energy diagrams contributing to order i^. Wiggly/dotted/full lines denote 
gluons/ghosts/scalars, respectively. In the notation of Eq. (i27|l . the first four diagrams give 
n^^^'°, while the last is n^°^'\ 

5.2 Two-loop gap equation 

At order £'^, the gap equation Eq. (p6]) reads 

= u''^\-m^)(l + dp2U''^\p^)\p2=_^2)+ui}\-m^) (32) 

= (<)'^ + n¥)'°) (i + a,.n¥)'° + 5,.<)'^) + <)-° + n«'^ 

where in the last line we have used Eq. (j30p and grouped together terms which will prove 
to be gauge invariant. There are 38 genuine two-loop diagrams contributing to 11^ , shown 
in Fig. [2j These can be expressed in terms of six scalar master integrals. Note that in 
unitary gauge (^ — )• c«, to be taken before regularisation) the ghosts and pseudo-goldstones 
decouple leaving only nine diagrams (see also |16j). The one- loop diagrams with one counter 
term insertion are shown in Fig. [3j A tree-level diagram with two counter term insertions 
is not one-particle-irreducible and hence does not contribute. From appendix O (where the 
renormalised 3d coupling ^^(/i) = /^"^'^S'bare ^^^ introduced), the different contributions to 
the two-loop gap equation are 



r(2),0 , ^(i),0, „^(i),0 ^^^2e fg^\ /^ _ ^ ^ ^^ ^2 

Svr / V20e 10 10 Am^ 



iif>- + nra,^nr "-™ (%ii) (^-io.6452 + i + ^i„£,) +oW, (33) 



T Svr VS 2 8 2V?-1 2^^y ^'' 

^ ^ 8^m V 32 8 8 2^ - 1 2^^) ^ ' ' 

where all quantities are understood on-shell (p^ = —m?) and for an analytic expression of the 
two- loop constant 10.6452 we refer to Eq. ()62p . Note that in general the on-shell self-energy is 
gauge dependent. However, the parts of the gap equation pertaining to the non-linear sigma 
model without counter terms, i.e. the first of Eqs. (133p . as well as the sum of all counter term 
contributions are separately gauge invariant, thus leading to a gauge invariant solution for 
the pole mass in the resummed theory. Our result for the second line of Eq. (j33p corrects an 
error in an earlier calculation for SU(2) [16j, which led to a ^-dependent pole. 
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2 '^ S/* 2 '^ }v* 'v{ )n/» 2 '^ 5^^ 2 
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Figure 2: The 38 diagrams contributing to U^ . Notation as in Fig. [TJ 

Figure 3: The 7 diagrams contributing to Ilrp'' . Notation as in Fig. [TJ 

The first of Eqs. p3p features a divergence as d — )■ 3, which is removed by mass renor- 
mahsation according to rn^^j.^ = m?{fj,) + 6m?, where Sm^ = — ■^— ( ^ ^ ) with x = 1 



20e \ Stt 
,2/ 



(x = 47re '^) specifying the MS (MS) scheme, respectively, and m (n) denotes the renor- 
mahsed mass. The renormahsed two-loop gap equation reads (with fl^ = 47re~'^^^) 

= m^ + —^ o 7-T7Tln3 +M^ -10.6452 + — + — In -^ . (34) 



Svr 2 V4 16 / V 8vr y V 10 10 Axm , 

Note that the renormalisation prescription has introduced scheme as well as scale dependence. 
In the following we pick the MS scheme (x = 1). The scale dependence is formally of 
higher order: since the bare mass is scale independent, the renormahsed mass reacts to a 

scale variation as m?{fi) = m^(^o) ~ f ( %^ ) ^^{f^/f^o)- ^^ ^ truncated perturbative series, 
however, this scale dependence remains and can be taken as an estimate for the size of higher 



q^N 


0.1 


0.3 


1 


3 


10 


Ci 


0.1692 


0.1651 


0.1605 


0.1562 


0.1512 


C2 


4.4 • 10-9 


1.3 ■ 10-^ 


4.4 • 10-8 


1.3-10-^ 


4.4 • 10-^ 



Table 1: Scale-dependent solutions of the two-loop gap equation Eq. (p4|) . m = Cig N 



order corrections. For a particular choice of scale, fi = m, the logarithm can be absorbed into 
the pole mass, such that the gap equation reduces to = rn^ — 0.07114g'^ Nm — 0.015165^A^^, 
with positive solution m2-ioop = 0.1637 5^ A''. 

Eq. (|34p possesses non-trivial, real solutions which are listed in Table [T] (again in the MS 
scheme, x = 1; the second solution, C2, is almost always close to zero, and we shall hence 
still call it trivial). Since the gluon does not correspond to an asymptotic particle state, the 
scale dependence of its pole mass is expected and in complete analogy to the two-loop pole 
masses of the electroweak gauge bosons [21]. Note that the pole mass changes by < 10% only 
as the renormalisation scale is varied over two orders of magnitude. Together with the fact 
that the two-loop contribution constitutes a ~ 15% correction to the leading-order one-loop 
result, this points to a reasonable convergence of the resummation scheme. 



6 Conclusions 



We have generalised a non-perturbative resummation scheme for three-dimensional Yang- 
Mills theory based on the non-linear sigma model [8| to SU(A^). Adding and subtracting 
a covariantly coupled scalar field allows for a gauge invariant gluon mass term regulating 
infrared divergences encountered in bare perturbation theory. We have established that this 
leads to gauge invariant physical results by analysing the BRS-invariance of the resummed 
theory. As an application, we have calculated the transverse gluon propagator and evaluated 
its pole mass by means of a gap equation, which we explicitly verified to be gauge invariant 
through two loops, thus correcting an error in [16J. The pole mass requires normalisation. We 
have employed the minimal subtraction (MS) scheme, through which it acquires a weak scale 
dependence. We find the two-loop correction to be ~ 15% of the leading one-loop result, 
and the scale dependence ~ 10% of the two-loop result when the renormalisation scale is 
varied over two orders of magnitude. Together, these two features might be indicative of a 
reasonable convergence behaviour of our resummation scheme. 

As a further application, the scheme lends itself to an evaluation of the (/^-contribution to 
the thermodynamic pressure in four-dimensional gauge theories. Preliminary investigations 
up to two loops have been reported in [22], a three- loop calculation is currently under way. 
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A Details of the calculation 



Due to the somewhat non-standard action with potentially high-order vertices, we automat- 
ically generate Feynman rules as well as a model file directly from Eq. (jlip . For two- loop 
self-energies, we potentially need vertices with up to six legs. 

For notational simplicity, let us write the loop expansion of the bare on-shell transverse 
self-energy as defined in Eq. ()23p as well as its derivatives as 



,2-^2n^^-(/) 



dp^ dm 



[m 



2\l-a-b \r-^ 
n>l 



"bare 



NJ{d,m) 



771^(1 — d) 



n 



(n) 



(35) 



where the ft'-"') are dimensionless functions of d, ^ and N only which are computed from 
n-loop Feynman diagrams, and J{d, m) is the massive one-loop tadpole integral defined by 
Eq. (|64p . Note that the square bracket in dimensional regularisation expands as 



"bare 



NJ{d,m) 



m?{l — d) 



=3-2e 



Sirm \2mJ 



2e 



[l + 3e + 0{e^)) , 



(36) 



fi ^"^^barc ^^"^ ^^^ usual MS scalc 



with renormalised coupling (note that Z„2 = 1 in 3d) g 

p? = 4iTe~^fi'^. From Eq. (j35p we explicitly have, in the notation of the main texlo, 



n, 



(i),o 



n, 



(2),0 



(-m^) 



(-m^) 



m 



m 



ffbare ^ JJ^, m) 

m,2(l — d) 
NJ{d,m) 



n, 



(1) 

00 ' 



"bare 



dp2llL {p 



Tm?(l — d) 

gl,^^NJ{d,m) 
m?{l — d) 



fT(2) 



n 



00 ' 



n 



(1) 



10 ' 



n, 



(i),i 



(-m^) 



-m 



"bare 



NJ{d,m) 



771^(1 — d) 



n 



(1) 

01 



(37) 
(38) 
(39) 
(40) 



From here, the calculation proceeds via standard methods. All diagrams that we need are 
generated with QGRAF [23J. We then shift momenta to our conventions, apply colour and 
Lorentz projectors, perform colour traces via the Fierz-identity, rewrite scalar products in 
terms of inverse propagators, and perform derivatives for Ilab on the integrand level. Using 



The last of the four relations is non-trivial and follows from realising the mass-shift m — >■ {l~£)m needed 
for the resummed theory Eq. (|11[) by a translation operator exp{p^£d^2) followed by the on-shell condition. 



10 



^, we obtain an intermediate result for the H in terms 



finally the on-shell condition p^ = —m 

of dimensionless one- and two-loop on-shell integrals I, defined as 

i{si,...,S4) = 
/(si,...,sio) = 



1 r d'^k I 1 




J{d,l) J (2^)'^ [k^ + ssY^ [{k-pf + s^Y- 

1 r d'^ki r d'^k2 1 1 


p2 = - 


[J(d, 1)]' J (2vr)^ J (27r)'^ [kf + s^Y^ [k^ + 
1 1 


siY^ 
1 



(41) 



.(42) 



[{ki - /C2)2 + SsY' [{kl - pf + S^Y^ [{k2 - pf + Sio]^^ 

where the normalisation factor J is a one-loop massive tadpole as defined in Eq. ^A 

In a next step, using symmetries and reduction relations following from systematic use 
[Ml [25] of integration-by-parts (IBP) identities [26], we arrive at d-dimensional expressions 
in terms of a few master integrals, as listed below. 

A.l One- loop computations 

Applying the necessary projectors on the sum of diagrams depicted in Fig. [1] as well as its 
p^-derivative (which we take at the integrand level) and performing a reduction to master 
integrals, we obtain the d-dimensional results 



(1) 

00 

(1) 

10 



ai Kl +a2K2 , 

bi Kl + b2K2 + bs K'l + 64 i^2 



(43) 
(44) 



with master integrals Ki listed in App.[B]and coefficients 



9 



Ol 



{4d - 5) 



1 



bi = -{d-2){Ad-5) 
lb 

^(1-40, 



02 
62 



{2d - 3){2d - 5) 



4 

-(12^2 _3i(i + i8) , 

8 

1 



(45) 
(46) 

63 = T (1 - 40 , &4 = - (3 - 2d) . (47) 

Similarly, either from the diagrams of Fig. [3] and using Eq. (j40p in reverse, or directly from 
the m^-derivative (easily taken at the integrand level) of the first four diagrams of Fig. [H 

4biKi + b5K2 + b3K'i + b4K'2, (48) 

65 = - {d^ - M^ + 4d - 3) . (49) 



n(i) 



Note that IIqq is gauge parameter independent, while II\q and IIqi are not. However, their 



difference 11 



n 



(1) 
01 



-^— IIqq is gauge invariant. Let us remark that this relation between 



the three II'^^ is not just a coincidence, but the one-loop on-shell case of a general relation 
which follows from using that dimensional analysis gives IIt ~ mass^ and g^ ~ mass^"'^, such 
that 

p^dp2+m^d^2 + ^^g^dg2] UT{p^,m^,g^) = l-UT{p^m\g^) . (50) 
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A. 2 Two-loop computations 

Applying the necessary projectors on the sum of diagrams depicted in Fig. [2] and performing 
a reduction to master integrals, we obtain the d-diniensional result 

n{,J = C1K3 + C2 K4 + C3K5 + a Ke + C5 Ki Ki + C6 Ki K2 + c^ K2 K2 

- ai 63 Ki K[ - 02 63 K2 K[ - ai 64 Ki K'2 - 02 64 K2 K'2 , (51) 

with master integrals Ki listed in App.[B]and coefficients 

ci = ^(d-l)(176d-245), (52) 

C2 = -— {U4d^ -712cf + 1241(1-760) , (53) 

64 

lOSOOd^ - 70632ff3 + I65227(i2 - 166654d + 61752 . ^, 

^^ = 192(3d - 4) ' (^^) 

C4 = -(d-2)(32ci3 -312(i2 + 656(i-405) , (55) 

64 

C5 = T^ (32d^ - U8d + 155) , (56) 

3 

ce = (16(i^ - 188d^ + 668(i2 - 940(i + 465) , (57) 

16 

C7 = -— ^^^ (24^5- 164d^ + 452^3 -680(i2 + 597d- 242) . (58) 

Individual diagrams do have contributions to IIqq that are proportional to 1/A^^ or 1/N'^, 
but these cancel in the sum, leaving IIqq A^-independent. 

Note that IIqq + IIqq IIj^q as well as IIqq + IIqq IIq]^ are gauge parameter independent. 

A. 3 Results in 3d 



Let us here collect the expansions around d = 3 — 2e of Eqs. (|43l) . (I44j) and (1481) 

n™ '■"' _5^1n3 + 5 + 0(£) =s -3.57579 + 0(£) (59) 

1° 32 8 8 2^/C - 1 2 ^ ^ ^ ^ ^ ^ 

°i 8 2 8 2^/^ -1 2^^ ^ ' ^ ^ 

satisfying IIj^q — IIq-^ ss — ^IIqq + 0{e), as well as of Eq. ([5T|) (subtracting ^-dependence) 

ft£'+n«<' '*" A + ?f2/(lZ!L_M?)_l^f_^ + 6U,(l/3)-2L,,(l/9) 

00 00 10 20e 32 ^ 512 V 6 ^ ' ^ ^ ' ' 

17069 16761 „ 369, „ 9 , 2o .^/ n 

+ + In 2 In 3 In^ 3 + O e (62) 

4800 320 8 512 ^ ^ ^ ' 

« ^- 10.6452 + 0(e). (63) 
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frW 



(1) 



For comparison, U\}; = 47r/2(0 + 0{e) in the notation of Eq. (21) in [16], but H^/ / 
— 47r/i(^) + 0{e), pointing to an error in that reference. 



B Master integrals 

We hke to work with dimensionless and measure independent integrals, so let us normalise 
each loop by the massive one-loop tadpole integral J{d,m), which for our choice of measure 
reads 



J{d,m) = J 



d'^k 



1 



l,("^Y\(r-,/2) 



d=3-2e m f ne^'^ \ , .^, ')-.\ 

'l + 2e + 0{e^)) . 



4tt V rn-2 



(64) 



Prom reference [2Yj (and using Eq. (1) therein as well as Eqs. (1,3,9-11,14) of [28j for K^), 
we get expansions around d = 3 — 2e for all on-shell master integrals that we need: 



Kt 



K[ 



K'l 



K:, 



m 



0/4 



1(1,1,1,1) 



d=3-2e In 3 



+ 0{e), 



^'-^/^J^^ = ^ai>e,e) 



.^_^2. li^2v|+l+0.. 



2 2V?-1 



m 



-Q^/- 



j =1(1,1,1,0 



d=3-2. 1 Vg + 2 

- --ln—^ + 0{e), 



m 



<I>/^= 



d=^2e /(1/3) - /(7/9) 
V2 



= 1(1,1,1,1,1,1,1,1,1,1) 

OS 

+ 0{€) with /(x) = 9Li2(x-F Vl-x2) 



167r2 X 0.000245310499 • • • -F ©(e) , 



i^4 



[^4 _Q_ / j2j ^^ ^ ^(^^ -^^ ^^ Q^ ^^ -^^ ^^ -^^ 0^ ^^ 



d=3-2e 1 



vr 



In^ 3 - — + 6Li2(l/3) - 2Li2(l/9) 
6 



+ Oie) 



Wtt^ X 0.00121156- •• +0(e) 



K. 



m 



^/j' 



7(1,0,1,0,1,1,0,1,0,1) 



d=3-2e 1 /I , , ,^, , 

^ 4^+(2-21^2)+0(6). 



(65) 
(66) 

(67) 

(68) 
(69) 
(70) 
(71) 

(72) 

(73) 
(74) 

(75) 

(76) 
(77) 

(78) 



13 



Two trivial one-loop massive vacuum master integrals read 

K2 = \ (~) /j] =1(1,0,1,0) = 1, (79) 

K'2 = [_Q_/J] =/(l,0,C,0) = e('^-2)/2. (80) 

One non-trivial fully massive vacuum master integral, expanded around d = 3 — 2e, reads 
(cf. reference |29j): 

Ke = [m2Q/j2] =7(1,1,1,0,0,1,1,1,0,0) (81) 

d=3-2e 1/1 3\ ,^, , 1 ,^, , , , 

« —- (- + ln-j +0(e) «— -0.9054651081 +0(e) . (82) 
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